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ABSTRACT 


A  finite  eleaent  model  is  developed  for  the  Class  IV  flextensional 
underwater  transducer  shell.  The  stiffness  and  mass  matrices  for  a 
finite  cylindrical  type  shell  element  of  linearly  varying  radius  of 
curvature  are  presented.  The  element  has  36  degrees  of  freedom 
corresponding  to  the  nine  generalized  coordinates  1/^ 

at  each  node.  All  six  rigid  boJy  modes  for  the  element 
are  adequately  represented  by  this  model.  The  element  is  used  to 
predict  the  free  vibrations  for  the  Class  IV  flextensional  underwater 
transducer  shell,  and  the  results  are  verified  experimentally. 


IMTROIXJCTION 


The  purpose  of  this  investigation  is  to  develop  a  finite  element 
model  for  the  Class  IV  flextensional  underwater  transducer  shell.  A 
picture  of  a  Class  IV  transducer  is  shown  in  Figure  1.  The  Class  IV 
flextensional  transducer  design  was  originally  proposed  as  having 
possible  applications  as  a  sonobuoy  transducer.  A  finite  difference 
model  of  the  Class  IV  transducer  has  been  derived  by  Rutledge  and 
Royster  (l),  (2),  and  (3).  Due  to  the  fact  that  plate  effects  were 
neglected  by  Rutledge  and  Royster,  an  attempt  is  made  to  develop  a  more 
complete  and  accurate  analytical  model.  This  paper  describes  such  an 
improved  linear  analytical  model  developed  for  the  calculation  of 
axisynmetric  modes  of  vibration  and  natural  frequencies  of  a  cylindrical 
type  (oval)  shell. 

A  finite  element  technique  is  utilized  to  construct  the  total  shell 
stiffness  and  mass  matrix  by  subdividing  the  shell  into  a  set  of 
axisynmetric  shell  components.  Novozhilov*s  "Thin  Shell  Theory”  (4) 
is  utilized  for  describing  the  state  of  strains  and  stresses.  The 
stiffness  and  mass  matrices  for  the  shell  are  obtained  by  superposition 
of  the  stiffness  and  mass  matrices  of  the  individual  shell  element  (5) 
which  are  computed  by  using  a  Ritz  approach.  The  superposition  technique 
assures  displacement  compatibility  and  force  equilibrium  at  the  interface 
between  components.  After  the  systems  for  stiffness  and  mass  matrices 
have  been  formulated,  displacement  boundary  conditions  are  introduced 
by  removing  appropriate  rows  and  columns  corresponding  to  points  on  the 
shell  which  are  rigidly  restrained  from  motion.  The  shell. 's  natural 
frequencies  and  node  shapes  are  obtained  from  the  eigenvalue  equation 
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constructed  with  the  total  stiffness  and  mass  matrices . 

in  which  GJ  is  the  circular  frequency  of  the  system  and  is 
the  modal  vector  where  the  nine  generalized  coordinate  component  of 
each  node  for  all  nodal  points  are  represented. 
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Am^YSlCM  MODEL 


Displacement  function 

As  pointed  out  on  page  ^2  of  reference  (5)»  the  first  requironent  of 
a  finite  element  is  that  it  must  he  capable  of  undergoing  rigid  body 
motions  with  little  or  no  strain.  Hence,  the  cylindrical  type  shell 
element  will  require  special  displacement  functions  in  order  to  ensure 
that  all  six  rigid  body  modes  are  included.  First,  it  is  clear  that  the 
twelve  term  polynomial  in  x  smd  y  used  for  the  plate  bending  element  may 
be  used  here  for  the  radial,  circumferential,  and  meridional  displacements 
WjV.u  (figure  2). 

'•■ft:)) .. -.Only 

«ii  x'+ OX) 

+0,^x7*; 


The  36  constants  in  equation  (l)  to  (3)  are  determined  as  functions  of 
the  36  modal  displacements  Ui,>rc^Wc^  aji, 

for  the  element  ^  *  "^^se  are  then  substituted  back 

in  equations  (l)  to  (3)  yielding  expressions  which  completely  difine  the 
displacement  anywhere  on  the  element  in  terms  of  the  modal  displacements. 


'X*-  C 


Figure  'd. 


Finl  •  jyli  i; 

;?ric  i  .^e  shell  element  of  linearly  varying 
of  curvature 

Point  1 

,  Y=.o 

Point  2 

x«a  ,  ysO 

iVWxlVi 

Point  3 

xra  , 

'*.*'<«»  .Mliail®).) 

Point  4 

X=0  ;yrb 

Therefore;  the  36  constants  are  determined  easily  from  the  above 
stated  36  equations. 
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■  -is. 

'^\''-'.  'N  '  \ 

''•i'  ^■^■'ilj^'^syinme'try  \if  Cl^  where  i^v412^is  given  by  a  function 

;  '‘•'^■i^^iTSi.-Jn,-aiXi,Uy,.,H^,Uxi,  >y.,\L^,J.^^,\J. 


V .  X 


V4 


N 

•I^J 


WIIX"' 

and  . 


vJv^A)y■v■^^l|^>^5 

The  element  with  4  nodes  1.2,3.  and  4  thus  possesses  36  degrees 


of  freedom,  determined  by  the  element  displacements. 

% 


ere. 


£ 


et. 


Ni 


The  displacement  within  the  element  is  uniquel.;^’  determined  by  the 

isl" 


nodal  displacements 


and  the  position  x  and  y. 


where  matrix  is  gwen  in  Appendix  A. 
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Strains  8md  Stresses 

From  Novoshilov’s  "Thin  Shell  Theory"  (4),  the  strains  and  stresses  for 
the  thin  orthotropic  shell  and  for  the  case  of  an  arbitrary  radius  of 
curvature  shell  is  as  follows 
Strains 

'  0»3- 

^  - 1.  The  ,  1C 

'  Ox*  / 

Stresses 

*■  C-£)c+VCy) 

S'  i!f«y 

A  general  elastic  behavior  is  assumed,  and  the  relationship  between 
stresses  and  strains  is  linear  and  of  the  form 


Mil 


56X 

) 


where  Ci>]  is  the  eleistic  matrix  containing  the  appropriate  material 


properties  and  is  given  in  Appendix  A. 


Since  matrix  M  is  already  calculated, 

then 

Matrix  is  easily  obtained  and  is  given  in  Appendix  A. 


The  stiffness  matrix  becomes 


-1^  o  o 

and  the  mass  matrix  becomes 


b  > 


X  ^ 

""72  O  0 

where  P  is  the  density  of  the  shell. 

a 

A  5  point  Gauss -legendre  numerical  integration  is  carried  out 
through  the  use  of  an  IBM  370/l65  computer  (6).  The  stiffness  and 
mass  matrices  of  the  element  will  clearly  always  be  square  and  of 
the  form 


Kk 

"111  "'ll.  MbMiy' 

^2l  ^2  1^2+ 

r.  -1^ 

"•«  "to.«hsMn. 

^31.  K^5  \</^4 

•V. 

14j  Ktx 

%  *»»n.wv(s 
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in  which  k^,  .  ,  •  etc.  and  •  etc.  are  sabmatrices  which  are  again 

square  and  of  the  size  9x9.  Then  the  totial  stiffness  LK] and  total 

CM]  Of  the  shell  are  obi;ained  by  proper  superpoistion  of  the  above 
mentioned  submatrices.  That  is 


The  "inside"  summation  is  now  taken  over  all  the  elements  of  the  shell. 
If  a  particuiar  element  does  not  in  fact  include  the  node  in  question, 
it  will  contain  no  submtvices  with  an  i  suffix,  and  therefore, 
its  contribution  will  simply  be  zero. 


After  introducing  the  boundary  conditions  to  both  matrices  Ck]  and 
by  removing  the  appropriate  columns  and  rows  from  them,  then  the 
problem  is  reduced  to  the  simple  form 


If  the  [K3  and  Q^3  matrices  are  of  order  nxH  ,  then  different  values  of  the 
angular  frequency  fij  can  be  computed,  and  each  of  these  frequencies 


will  be  associated  with  a  particular  node  vector  which  also 
can  be  computed  by  utilizing  a  digital  computer  (6). 


DKCUSSION 


The  approach  adopted  here  is  known  as  the  displacement  formulation 
which  is  equivalent  to  the  minimi7,ation  of  the  total  potential  energy  of 
the  system  in  terms  of  a  prescrioed  displacement  field.  If  this  displace¬ 
ment  field  is  defined  in  a  suitable  way,  then  convergence  to  the  correct, 
result  must  occur.  The  process  is  the  equivalent  tc  the  Ritz  procedure, 

8ind  as  a  consequence,  the  so  obtained  approximate  solution  is  eJ-ways  an 
upper  bound  to  the  exact  solution  which  is  very  much  desired.  It 
should  be  emphasized  that  this  bo'jr.d  is  not  available  for  a  solutions  in 
which  discontinuity  at  interfaces  arises  (5). 

The  major  characteristics  of  the  presented  analytical  model  formulation 
are:  l)  an  accurate  representation  of  the  shell  geometry,  2)  the 
dynamic  mass  matrix  representation,  and  3)  the  requirement  of  rigid 
body  displacements. 

The  model  described  herein  has  been  utilized  to  calculate  eigenvalues 
and  eigenvectors  of  the  Class  IV  flextensional  underwater  transducer  shell 
(6).  Using  a  mess  of  6  elements  as  described  in  reference  (6),  the  first 
natural  frequency  (1295.25  cycles/sec.)  is  calculated  within  3*62^  as  compared 
to  the  experimental  value  (1250  cycles/sec).  Using  a  mess  of  8  elements  (hK2) 
the  first  natural  frequency  (1289,73  cycles/sec.)  is  calculated  within  3.17^ 
from  the  experimental  value.  This  clearly  indicates  the  accuracy  and  rapid 
convergence  of  the  herein  developed  finite  element  model. 
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N  (1,Q2)  =  N3  tA^3)=M GA-O-  ^  (r'^'^  t'X^Y) 

N  (1^0  -  lY  Q-,^^>)^CV7>  b  C-^y^  X  y!) 

N  Uj,a.8 )  =  iJ  LA-A‘=')-tY  C>;3o)  =•  xy  S'xV  -  b'XY^ 

cSlY*-(-  ^xN  4  S'XY^ 

N  ^  cp^ba-y  =  M  c^A'syr-  a.(xy  _  a,x?Y  tx^Y) 

=  -  WC^G6)=  b(_y^  y.'i'^f-y.'i^) 

'wWrt  -AjL  o(  -Hv€.  UC^TTi  f'T'O-V  ^  \  0€ 

o.\aoY^  cliX  ^^ucJi  4-0  ^ero  , 
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ns-x-^y^Y-a-y^) ; 

'Tba,'2.)  '^Ci:i')T 

BKf-)  -T^ i^.(i-XY-3'/<-3,Y'^Bx\i- iy,f'tS.'<\tL'^- 

SX^Y-SX^^)  i  •3C3,a)ri<l+^)(-Y-!iX+e.XY+3)X?^!iX- 

:5Ci,5)=  -(^.(-i-vY  +  a^-axY)^  &(jl;5)^^('l- 

X  Y  -  3X^3^'"^+  5X*Y  +■  B)/X%  a,x^  5 Y*- sx^'  -  axY^3 - 

^  ■(-i+X+a.Y-a.XY)  J  -803,3)  =  -^.(-i+iiXTGV-C.xW) 

f)  -  i  -Y-  4X+  -l-x Y  +3Yi^3'XV  J  3Ci^4)  :r  O 

'BC3,+.)=  ^•C-X^l.X-TC^), 

m,r)  =0^3. 04.  S’)  ( 

x-xY-ax’^+sx^ftx^-x^v) ;  @c^s-)  r  (i+*?|)(i-)-fx+ 

4XY  +  3X^-3X^y)^ 

■?)CL>)--±_  (-4-T4Y+iiX-C=.xY)  ;  3C4,fc)  -  2r,(_x-XY  -S^u 
SxY  +  X^xN;  ^  6C3,t)  =  -^.(;L-4-X-3X^), 

,  '^0V)=o, 

X  -4Y  i-4-xY  -r3V^-3XY'^ 
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'BCi,q)«o  ,  3C'5-,^)=4.(v-xY-aY-^xY^y!>y^3- 

f-(-++f-X-^feY'-  £XY)  j  BCa'^)3.-|:t(^.t4^-Y-3y^3^ 
S(l/0)  =  ^CYt-S>X-£.XY-3V-4xS4xN+aY’^)^QC2;,o):.0^ 

■&C2-,l'2.)= 

'  '^~  4’  ^  ^  -3xy  -  a  x^-t-axV  ■+  a.xY^ ) , 

■  ^  0  -  ^-  C^  ■*•  ^)  CY^£.x-&.xy-3,Y-£.x^  £y V  ^ 

Bci,a)  =-Gi  (i-v-A/+2xv)^  6c2^ia)=l(xY-T3y^- 
-  3xyl  ay2^ax2^  +axY")  (-x-^axy) 

'feC.i.'l^y  ^  -SX+Syy-r'iX^-’ixN^ 

BC-f-,  i  +)  -  O  ,  B cx,  i4) t:  ^  (Ivl ) (_xV)  -  C 

-X  +  x\  T  X  -  xN")  '£>C3,  1 4)=  (uai.)(_a-x-4-axY  4.BX  ~ 

3xN), 

■Bci,-(nr  -  Sii.(_i4Y43,x-s,xY:)  ^  ac^,i  s')^^(-v^x"‘y»xVv) 


C1  -S'/ +V^)  ,  3Ci,tC0  -Oj  'B©/b>X-f)(V+3XY'‘ 

'S>a,i7)  r^o_,-s>  cyi)  n  (u^X>^-4xv+3>>y  -■^^(x'y-axv> 
Xf)^55C^,'7)= 

'Bcij  &)  ^  83  .  A(xY-gXY^xV^)  _  ^(-4-x*&x7)  ^ 

■fees, 183= -^Ql-4-Yf8'y’-3 

j  8Ca-,H3  j 

Bca,  I3)  =  ^  C-XT37^%&xy-  5xIgxY^)  , 

'^Cl,ao3  =  o,  "Bca^^o^  ^(uD(-X4.3X^S,yY-?lx-(;X7^3  -- 

•  "Y+&XY +3Y*^(oxVSyO; 

'Bax0  =  -^^C'^-Sxx3  ,'K.a,^lO-^(-XY+3xV3'XY'^ax^Y- 

cC*"  ^ 

-  %C^-a,xV)^3Ci,aD=  -S:|(-i4!ix+<iY-tx":i./) 

D  CL  o 

'?ci,>^>?)=  -axYt^.x^'  j  B a,aa) = b<8,s?3.3  -eC-x^ ^ 
BCi  a.3)=o^  Ba/13)=  ^(i+l-U-xVi-^^^x^TxN)^ 
BC3,530=  Ci^^)C-^V+3x"v)^ 

•Ba,a4Y^-^(3aY+fcxv)  j  sa,3^=^(^x^4.^Y), 

BC3,a+)  ■=  -  ^(-a>;-3x\ 
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,  60^,5SO  =  0,  BC^,,aO=-axV^3>Y^  - 

a7>o,  «c>^a7>k  (-x7rX/J-|c-av?6x  b(?,s-7}- 

-|p(^-av^3v^3^ 

'3Ci,a8;)=:  i(7-i,xv  -2>7+&,<V  ’•avOj  aca^as) =q,  sci^a.8i= 
^(.^+  (S  y-3X^  c=xv  -  G'i^t^yX &y  7^)^ 

'Ba,aot:)=(3^  £>Gi:,'^=^:>=^Ci+4)Cx+&7-a'<-&xv-c.v^ 

4  X  V^_)  -  '^^(ri  -Psi-iy^i- 3  X  ) , 

BC3,a^3  =  ic  i  )C  V-  4XV.-av\  fc  V^V  X  A  Y^X 

.  '^C-l,l>oJ)  =  -  ^  (_- "(  +Uxv  )  ^  6a, bo)  =  ^Cx'7  +37 -3>xX  - 
.  3XY  — ^y'+Sx^Y  +S.aY"*  J  .3LB,Boj  r a-4X’*!<jY'*‘4x +47^) 

^  CLb 

'BCi,'^0  -  ^-4-xY  +3xV  "aa,  3.^-) -  o_^  Bc=i,B.'a.)-“(x-ax^x^^ 

'BC1-,32)=  O  ^TB-Ci/sZj- &(a£Xx-AX+X^)-'^'"^CXY- 

SixX+xXx  e(Ba)=(i+|sxW-Yxy+3x^YV  ^ 
Xi33)^-|(-tY  +  4xv')  ^  £,t23)+^C^V-?^xVxS)^ 
BC3,BB)=-^(i-Y-X-3xX 

b 

'feCMtj  *  '^ci'i')  J  T-a;iA-)^o^  2o^3+>  _aY+axv+. 
'5v\.3xv^ 

•  V 


'&ca>r:)=  Ci+^X-^v^a.)<Y*^v-3xVt) 

^■^(-Y +xy^y"-xv^).  ^Ci+^)C 

Y-y"3, 

3.(_i^3<i; =q^  B.ca,3!y  ^  (-  T  V -X  Y^j)  - ^ 

2Xv'^i^V-GXv3  j  "" ^  # 


